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1.  INTRODUCTION 


In  this  report  we  summarize  the  research  effort,  funded  under  ONR  Grant 
N00014-91-J-1004,  on  probabilistic  modeling  and  statistical  inference  for  random  fields 
and  space-time  processes.  The  research  we  have  pursued  in  this  project  consists  of  the 
investigation  of  a  set  of  interrelated  topics  involving  the  development  of  new 
mathematical  methods  for  the  challenging  problems  of  random  field  analysis  and 
inference  and  the  application  of  these  methods  to  problems  of  practical  significance.  In 
particular,  the  problems  of  interest  in  this  work  were  motivated  by  ^d  directly  adless 
issues  that  are  central  both  to  the  challenging  large-scale  data  assimilation  and  estimation 
problems  arising  in  physical  oceanography  and  to  a  number  of  other  remote  sensing  and 
imaging  problems  of  direct  interest  to  the  Navy. 

During  the  grant  period  we  have  had  what  we  feel  is  considerable  success  in  our 
research  efforts.  In  particular  the  research  supported  in  part  by  our  ONR  grant  has  led  to 
90  publications,  including  40  journal  papers,  36  conference  papers  and  presentations,  2 
book  chapters,  3  SM  theses,  and  9  Ph.D.  theses.  A  number  of  additional  papers  and  theses 
at  both  the  SM  and  Ph.D.  level  are  currently  in  progress.  Furthermore,  our  research  in 
this  area  has  received  considerable  national  and  international  recognition.  In  particular. 
Professor  Willsky  has  given  several  keynote  and  plenary  lectures  on  research  topics 
supported  by  our  ONR  grant.  These  include: 

•University  of  Southern  California  Signal  and  Image  Processing  Institute 
Distinguished  Lecture,  Feb.  1991. 

•Featured  invited  lecture  at  the  SIAM  Conference  on  Linear  Algebra  and  Its 
Applications,  Sept  1991. 

•Keynote  address  at  the  IEEE  International  Conference  on  Systems  Engineering, 
Aug.  1991. 

•Principal  lecture  on  Multiresolution  Methods  in  Image  Analysis  at  the  Tri- 
Service  Workshop  on  Statistical  Methods  in  Image  Processing,  May  1992. 

•Keynote  address  and  principal  lectures  at  the  inaugural  workshop  for  the  Centre 
for  Robust  and  Adaptive  Systems”,  Canberra,  Australia,  February,  1992, 
involving  both  researchers  and  members  of  the  Australian  government. 

•Featured  invited  lecture  on  challenges  in  signal  processing  at  the  SIAM 
Workshop  on  Systems  and  Control,  September  1992. 

•Keynote  address  on  challenges  in  signal  processing  at  the  Twenty-Fifth 
Anniversary  Celebration  for  INRIA  (Institut  National  de  Recherche  en 
Informatique  et  en  Automatique),  held  at  the  Ministry  of  Research,  Paris,  France, 
December  1992. 

•Plenary  address  at  the  1993  IEEE  Symposium  on  Image  and  Multidimensional 
Signal  Processing,  Cannes,  France,  September  1993. 

•Principal  invited  lecture  on  statistical  problems  in  image  analysis  and  random 
fields  at  The  Joint  Statistical  Meetings  of  The  American  Statistical  Association, 
the  Biometric  Society,  the  Statistical  Society  of  Canada,  and  the  Institute  of 
Mathematical  Statistics,  Toronto,  August  1994. 
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In  addition,  Mr.  Paul  Fieguth,  a  Ph.D.  student  with  Prof.  Willsky  has  had  several 
honors,  including  the  Student  Paper  Award  at  the  1994  IEEE  Oceans  '94  Meeting  in 
Brest,  France  and  invited  lectures  at  the  1994  International  Conference  on  Image 
Processing  in  Austin  Texas  and  at  the  upcoming  1995  International  Conference  on 
Acoustics,  Speech,  and  Signal  Processing. 

Furthermore,  during  this  time  period  we  have  had  considerable  interaction  with 
Navy  personnel  and  with  problems  of  direct  Navy  interest,  leading  both  to  the  application 
and  transitioning  of  our  research  and  to  the  formulation  and  study  of  a  variety  of  research 
problems  that  have  been  inspired  by  discussions  with  Navy  personnel  or  that  have  been 
motivated  by  Navy  missions  and  needs: 

•Our  initial  research  on  multiresolution  methods  for  estimating  velocity  fields  in 
space-time  data  and  sequential  imagery  has  been  applied  by  Alphatech,  Inc.  to  the 
problem  of  image-based  tracking  of  targets  under  an  SBIR  (Small  Business 
Innovative  Research)  Phase  I  project  sponsored  by  NAVAIR  and  directly 
monitored  and  supported  by  Dr.  Gary  Hewer  at  NWC  China  Lake.  The 
continuation  of  this  work  under  a  Phase  n  program  is  currently  under 
consideration,  with  the  endorsement  of  Capt.  Obar,  PMA  for  the  Sidewinder 
Program. 

•The  fi’amework  we  have  developed  for  modeling  random  fields  at  multiple 
resolutions  and  in  particular  for  estimation  and  detection  in  ffactal-like 
backgrounds  is  also  being  transitioned  by  Alphatech  in  the  context  of  detecting 
smaU  targets,  such  as  periscopes,  from  P-3  SAR  and  ISAR  imagery,  under  both  a 
second  NAVAIR  SBIR  and  directly  with  NWC  (Dr.  Oran  McNiel). 

•Prof,  Willsky  has  had  direct  contact  with  researchers  and  engineers  at  the  Naval 
Coastal  Systems  Station  (including  Drs.  E.  Moritz,  G.  Dobeck,  and  A.  Dubey) 
involved  in  the  development  of  multisensor,  multispectrd  methods  for  imaging 
and  automatic  detection  and  recognition  of  mines.  This  interaction  began  with 
discussions  at  the  1992  Tri-Service  Workshop  mentioned  previously  and  has 
continued  with  two  trips  that  Prof.  Willsky  has  taken  to  Panama  City  in  1993,  the 
first  at  the  invitation  of  Dr.  Dubey  (during  which  Prof.  Willsky  gave  a  lecture  on 
his  work  on  multiresolution  image  analysis)  and  the  second  at  the  ONR/CSS 
Workshop  held  in  November  1993. 

•In  addition.  Prof.  Willsky  has  interacted  with  researchers  at  The  Charles  Stark 
Draper  Laboratory  who  are  working  under  contract  to  CSS  in  this  area.  Indeed,  at 
this  time  Draper  and  Dr.  Dobeck  are  discussing  the  incorporation  of  some  of  our 
ideas  for  the  analysis  of  synthetic  aperture  sonar  data  collected  at  CSS. 

•The  algorithms  that  we  have  developed  for  multiresolution  sensor  data  analysis 
are  also  currently  being  transitioned  by  Draper  Lab  in  its  work  on  SAR  image 
analysis. 

•One  of  Prof.  Willsky’s  former  graduate  students.  Dr.  T.M.  Chin,  has  been 
working  with  Prof.  A.  Mariano  at  The  University  of  Miami  on  the  problem  of 
spatio-temporal  tracking  of  the  North  Wall  of  The  Gulf  Stream  using  methods  of 
optimal  estimation  developed  in  part  by  Dr.  Chin  during  his  graduate  work  at  MIT 
under  our  ONR  Grant.  The  problem  of  mapping  and  tracking  North  Wall  position 
is  of  considerable  interest  to  the  Navy  and  in  fact  such  maps  are  regularly 
produced  by  NRL  researchers  at  the  Stennis  facility  in  Mississippi. 
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•We  have  recently  demonstrated  the  power  of  our  multiresolution  random  field 
estimation  framework  in  the  context  of  processing  satelhte  altimetry  data  from  the 
TOPEX/Poseidon  satellite  for  the  estimation  of  ocean  height  perturbations  in  the 
North  Pacific.  Such  estimates  are  of  critical  importance  in  ocean  science 
programs  funded  by  the  Navy  and  by  other  national  scientific  organizations,  as 
they  provide  inputs  to  global  ocean  circulation  models  currently  under 
development  and  study. 

•Ms.  Lori  Belcastro,  another  of  Prof.  Willsky's  graduate  students,  is  the  recipient 
of  an  ONR  Graduate  Fellowship. 

A  complete  list  of  the  papers  arising  out  of  our  research  is  included  at  the  end  of 
this  report.  In  the  following  sections  we  provide  a  brief  summary  of  the  most  recent  of 
these  results. 

2,  MULTIRESOLUTION  MODELING,  ESTIMATION,  AND  STATISTICAL 
INFERENCE 

This  area,  which  represents  the  theoretical  basis  for  much  of  our  work  on  random 
fields,  covers  our  continuing  efforts  in  developing  and  exploiting  the  framework  we  have 
developed  for  multiresolution  modeling,  estimation,  and  statistical  inference  for  random 
fields.  Our  most  recent  publications  in  this  area  are  [58-60, 62, 71, 78,  80,  86].  In 
particular,  the  results  reported  in  [58]  provide  us  with  the  multiscale  statistical 
characterization  of  the  error  in  the  estimation  of  a  random  field  based  on  measurement 
data  that  may  be  sparse,  irregularly  sampled,  multiresolution,  etc.  The  significance  of 
this  result  is  twofold.  First,  it  allows  us  to  compute  error  statistics  for  random  field 
estimation  problems  of  considerable  size  (e.g.  1000  x  1000  random  fields),  something 
that  has  heretofore  been  impossible.  Secondly,  it  provides  us  with  one  of  the  key 
elements  needed  for  efficient  space-time  estimation  of  large-scale  random  fields,  since 
the  key  step  in  such  a  problem  is  the  use  of  new  data  to  estimate  the  error  in  the  estimate 
based  on  preceding  data  sets.  These  results  will  be  of  critical  importance  in  the 
continuation  of  our  research  on  space-time  data  assimilation. 

The  results  reported  in  [71, 78,  80, 86]  are  of  importance  as  they  provide  the 
statistical  basis  for  building  multiresolution  models  from  data  or  from  covariance 
estimates  of  a  random  field.  Given  the  substantial  advantages  of  our  multiresolution 
modeling  framework,  the  key  questions  are  in  identifying  under  what  conditions  it  is 
applicable-i.e.,  when  the  random  fields  of  interest  can  be  efficiently  modeled  within  our 
framework— and  in  developing  methods  to  construct  these  models.  Of  particular  interest 
here  are  the  results  in  [80]  and  [86].  In  particular,  the  method  described  in  [86]  provide 
an  extremely  powerful  methodology  for  building  multiresolution  models  Aat  have 
desirable  properties,  such  as  appropriate  levels  of  smoothness,  for  a  surprisingly  large 
class  of  random  fields  including  but  not  limited  to  those  that  are  typically  described  using 
Markov  Random  Fields  (MRFs).  The  work  in  [80]  demonstrate  the  promise  of  these 
modeling  methods  by  showing  their  applicability  in  modeling  S  AR  imagery  and  in  using 
the  resulting  models  in  discrimination  problems. 

3.  ESTIMATION  OF  GEOMETRIC  FEATURES  IN  RANDOM  FIELDS 

The  second  area  of  our  research  focuses  on  the  estimation  of  geometric  features  of 
random  fields.  Our  most  recent  work  in  this  area  is  described  in  [54-57, 61,  63,  64,  89]. 
The  research  described  in  [54-57, 61,  89]  deals  with  the  direct  estimation  of  geometric 
features  from  tomographic  measurement  data.  These  results,  and  in  particular  those  in 
[54],  provide  the  basis  for  a  new  approach  to  tomographic  reconstruction  from  extremely 
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limited  data,  as  is  typical  in  oceanacoustic  tomography.  In  particular,  our  results  indicate 
how  we  can  use  tomographic  data  to  estimate  the  spatial  moments  of  a  random  field 
which  in  turn  can  be  used  together  with  a  nonlinear,  divergence-based  regularization 
criterion  to  enhance  reconstructions  obtained  using  standard  tomographic  techniques. 
Thus,  for  example,  if  we  have  a  priori  information  about  the  location,  size,  or  shape  of 
spatial  features— such  as  cold-  or  warm-core  eddies— we  can  use  this  together  with  the 
estimated  moments  to  enhance  reconstruction  in  the  vicinity  of  the  eddy.  These  sarne 
methods  are  also  of  potential  value  for  nondestructive  evaluation  applications  in  which 
sparse  tomographic  data  are  to  be  used  to  detect  and  characterize  flaws,  cracks,  and 
corrosion.  The  work  in  [63, 64]  deals  with  the  problem  of  tracking  the  geometry  of 
objects  embedded  in  random  fields  that  change  with  time. 

4.  INVERSE  PROBLEMS 

The  third  component  of  our  research  deals  with  multiresolution  methods  for 
inverse  problems,  and  our  recent  work  in  this  area  is  described  in  [47-49, 65, 66,  68-70, 
72, 76-77,  81-84, 90]  In  particular,  the  research  in  [48-49,  65, 68, 72, 76,  84,  90]  deals 
with  the  development  of  multiresolution  methods  for  tomographic  reconstruction  that  (a) 
are  extremely  efficient;  (b)  allow  us  to  tradeoff  resolution  and  accuracy  in  a  statistically 
meaningful  way;  (c)  provide  direct  imaging  of  features  at  a  hierarchy  of  resolutions, 
allowing  us  to  identify  significant  anomalies;  and  (d)  provide  a  new  and  computationally 
superior  method  for  direct  texture  discrimination  fi-om  tomographic  data. 

Our  work  in  [47,  66, 69, 70, 77,  81-83]  deals  with  the  development  of 
multiresolution  inversion  methods  for  a  variety  of  inverse  problems,  including  those 
specified  by  the  partial  differential  equations  arising  in  inverse  scattering  problems. 
These  results  have  many  of  the  same  characteristics  as  our  results  on  tomographic 
reconstruction,  with  some  important  differences  in  emphasis  and  in  technique  due  to  the 
greater  complexity  of  the  Green's  functions  characterizing  inverse  scattering  problems. 

In  particular,  in  this  context  we  have  developed  a  precise  statistical  criterion  for 
determining  the  optimal  scale  of  reconstruction  at  each  point  in  an  image  and  for 
assessing  which  measurements  in  a  set  of  scattering  experinaents  provide  useful 
information  for  reconstruction  at  each  point  in  the  region  being  imaged.  In  this  way  we 
completely  expose  the  structure  of  multisensor  fusion  for  such  problems.  In  addition,  we 
have  now  also  developed  a  coarse-to-fine  approach  to  detecting  and  isolating  ^omalous 
regions  in  a  medium  from  inverse  scattering  measurements  [81],  allowing  us,  in  a  data- 
adaptive  way,  to  "zoom"  in  on  regions  of  interest  and  to  direct  the  limited  numbers  of 
degrees  of  freedom  in  our  reconstruction  into  areas  in  which  additional  detail  is  present 
and  can  be  estimated.  In  addition,  the  inverse  problems  of  inverse  scattering  are 
fundamentally  nonlinear,  as  are  many  inversion  problems  (including  oceanacoustic 
tomography).  In  much  of  our  work  in  this  area  we  focused  on  developing  inversion 
methods  based  on  linearization  techniques  (most  notably  the  so-caUed  Bom  linearization 
method), which  are  frequently  quite  accurate  (if  the  deviations  from  the  nominal 
background  represent  "weak"  scattering  perturbations).  In  our  most  recent  work  [83]  we 
have  begun  to  take  a  careful  look  at  the  fully  nonlinear  problem.  In  particular,  we  have 
developed  results  on  the  sensitivity  of  linearized  inversions  to  model  errors  induced  by 
nonlinear  effects.  One  of  the  novel  aspects  of  this  work  is  that  our  formalism  for 
determining  the  optimal  scale  for  reconstruction  allows  us  to  determine  the  r^uired 
resolution  to  which  we  must  accurately  capture  nonlinearities.  For  example,  if  the  scale 
to  which  statistically  significant  reconstruction  is  possible  is  coarser  than  the  scale  at 
which  nonlinear  effects  become  significant,  then  the  use  of  a  linearized  technique  is 
sufficient. 
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5.  APPLICATIONS  OF  MULTIRESOLUTION  STATISTICAL  METHODS  IN 
OCEANOGRAPHY  AND  COMPUTER  VISION 

Finally,  the  fourth  segment  of  our  recent  research  [67, 73-75, 79,  85-88]  deals 
with  the  application  of  our  multiresolution  methodology  for  estimating  random  fields  to 
several  problems  of  considerable  practical  importance.  In  particular,  the  estimation  of 
ocean  height  patterns  is  an  essential  element  in  understanding  ocean  circulation,  as 
variations  in  sea  level  provide  indications  of  pressure  differentials  that  dnve  the 
circulation.  Maps  of  sea  level  variations  are  often  used  as  an  end  result  in  themselves, 
providing  images  of  circulation  patterns— roughly  speaking  the  ocean  equivalent  of 
atmospheric  highs  and  lows— and  changes  in  those  patterns  in  time.  In  addition,  these 
maps  are  of  considerable  value  to  ocean  modelers.  In  particular,  global  circulation 
modelers,  such  as  our  colleagues  at  MIT,  use  such  sea  level  variations  both  to  drive  their 
models  and  to  correct  for  errors  between  the  predictions  provided  by  their  simulation 
models  and  the  estimates  of  ocean  height  based  on  measurements.  Roughly  speaking,  we 
can  think  of  this  as  an  extremely  large  time-recursive  distributed  parameter  estimation 
problem,  in  which  a  global  circulation  model  is  used  to  capture  ocean  dynamics  and  in 
which  maps  of  sea  level  variations  are  used  as  the  measurements.  In  order  for  this  to  be 
done  in  a  statistically  meaningful  way,  of  course  we  must  have  a  measure  of  the  quality 
of  such  maps.  That  is,  we  require  knowledge  of  the  covariance  of  the  errors  in  our 
estimate  of  ocean  height.  Such  covariance  estimates  are  of  use  for  other  reasons  as  well, 
as  they  allow  us  to  determine  if  particular  features  in  a  map  are  statistically  significant  or 
not,  a  point  on  which  we  will  comment  further  in  a  moment. 

There  are  a  number  of  significant  challenges  in  producing  such  maps  of  sea  level 
variations  and  their  associated  error  statistics.  The  first  is  the  fact  that  the  data  on  which 
such  estimates  are  to  be  based  are  sparse  and  irregularly  sampled,  s  is  certainly  the  case 
for  the  TOPEX/Poseidon  satellite  altimetry  data  that  we  have  used  in  our  work.  The 
irregularity  of  the  data  pattern  presents  a  major  challenge  as  there  is  no  regular  structure 
that  can  be  used  to  advantage,  thus,  for  example,  precluding  the  direct  use  of  efficient 
Fourier  transform  techniques.  Furthermore,  there  are  other  sources  of  nonstationarity  not 
only  in  the  expected  variability  in  sea  level  in  different  portions  of  the  ocean  but  ^so  in 
the  quality  of  the  measurements  provided  by  the  satellite.  In  particular,  the  satellite 
altimeter  provides  direct  measurements  of  the  distance  from  the  ocean  surface  to  the 
satellite,  while  what  is  desired  is  a  measurement  of  ocean  height  relative  to  the  geoid  (the 
equipotential  surface  of  the  earth's  gravity  field).  Since  the  geoid  is  not  known  perfectly, 
especially  when  one  accounts  for  the  presence  of  variations  in  ocean  floor  height  due  to 
underwater  trenches,  mountain  chains,  sea  mounts,  etc.,  errors  in  the  geoid  translate 
directly  into  errors  in  the  sea  level  estimates  produced  by  the  "corrected"  altinietry 
measurements.  Thus,  the  resulting  measurement  "noise'  in  the  altimetry  data  is  spatially 
nonstationary. 

In  addition  to  having  to  deal  with  irregular  data  and  nonstationary  phenomena,  we 
also  must  deal  with  the  fact  that  the  set  of  variables  to  be  estimated  is  extremely  large. 

For  example,  in  our  work  we  have  focused  on  estimation  in  the  North  Pacific,  a  region 
consisting  of  300,000  -  500,000  grid  cells  representing  the  finest  scale  at  which  we  may 
produce  estimates.  Thus  in  principle,  producing  a  map  of  sea  level  variations  over  this 
region  corresponds  to  a  300,000  -  500,000-dimensional  estimation  problem-and  on  top 
of  that  we  also  want  the  error  covariance  for  this  estimated  map!  To  our  knowledge  and 
that  of  our  oceanographer  colleagues,  no  methods  developed  by  others  are  capable  of 
dealing  with  problems  of  this  size.  In  contrast,  as  we  have  shown  in  [67 , 73-75, 79,  87], 
using  our  methods  we  can  perform  these  computations  with  considerable  ease,  allowing 
us  not  only  to  provide  the  maps  and  error  statistics  mentioned  previously  but  also  to 
investigate  a  variety  of  related  problems  of  considerable  importance  in  oceanography. 
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In  particular,  in  our  initial  studies  we  have  used  our  multiresolution  estimation 
methods  based  on  a  so-called  "fractal"  prior  model— i.e.,  a  l/f-like  spectrum— that  is 
consistent  both  with  the  statistics  of  the  satellite  data  and  with  the  form  of  the  models 
proposed  by  oceanographers.  Such  models  capture  spatial  structure  over  all  scales  and 
thus  are  well-matched  to  the  multiscale  nature  of  oceanic  phenomena.  In  addition,  as  we 
have  pointed  out,  the  accuracy  of  the  satellite  measurements  is  spatially  varying,  thanks 
to  the  varying  levels  of  uncertainty  in  the  geoid  correction  that  has  been  performed  on  the 
satellite  data.  We  have  captured  this  in  our  work  by  the  use  of  a  spatially- varying 
measurement  noise  variance  developed  in  consultation  with  our  oceanography  colleagues 
and  based  on  a  careful  examination  of  the  nature  of  geoid  errors. 

The  results  of  our  work  in  this  area  have  been  extremely  encouraging  and  have 
already  caught  the  attention  of  the  oceanography  community,  both  at  MIT  and 
internationally  (as  evidenced  by  Mr.  Fieguth's  paper  award  and  the  excellent  reception 
our  work  has  received).  In  particular,  our  method  allows  us  to  produce  maps  of  sea  level 
variations  at  far  finer  scales  than  those  that  oceanographers  had  previously  used.  In 
addition,  as  we  have  stated  previously,  one  of  the  major  challenges  in  statistical  inference 
and  probabilistic  data  assimilation  for  oceano^phic  analysis  is  the  enormous  size  of  the 
problems  that  must  be  considered.  Indeed,  using  standard  techniques,  such  as  Markov 
random  fields,  computing  the  estimates  for  the  ocean  region  considered  would  be  an 
exceedingly  difficult  computational  problem,  and  the  calculation  of  the  statistical 
characteristics  of  the  estimation  errors  would  be  prohibitively  complex.  Using  our 
methods,  the  processing  of  10  days’  worth  of  satellite  data  over  the  North  Pacific,  to 
produce  an  estimated  field  at  multiple  resolutions,  together  with  full  error  covariance 
information  takes  less  than  1  minute  on  a  Sparc  10  workstation. 

Thus  not  only  can  we  solve  problems  such  as  the  ones  that  we  have  just  described 
which  had  previously  been  considered  prohibitively  complex,  but  we  can  also  consider  a 
variety  of  deeper  and  more  ambitious  questions  of  considerable  interest  to  the  Navy  and 
the  oceanographic  community.  In  particular,  we  have  used  the  error  covariance 
calculations  to  allow  us  to  identify  areas  in  which  the  satellite  data  are  anomalous  and 
have  verified  that  indeed  these  outliers  correlate  remarkably  well  with  regions  in  which 
there  are  significant  bathymetric  features  and  geoid  uncertainty.  In  addition,  by  using  a 
higher-order  version  of  our  algorithm,  we  have  been  able  to  consider  the  joint  estimation 
of  sea  level  and  geoid  errors,  yielding  enhanced  maps  of  the  geoid  which  are  of  potential 
value  in  a  wide  variety  of  contexts. 

Moreover,  inspired  by  our  success  in  oceanography  and  by  the  fact  that 
oceanographers  are  actually  most  interested  in  the  gradient  of  the  ocean  surface  (as  this 
provides  information  on  currents  and  pressure  fields),  we  have  developed  new  general 
methods  for  surface  and  surface  gradient  reconstruction  from  noisy  and  possibly  sparse 
data  [85,  88].  These  results  are  not  only  of  interest  in  oceanography  but  also  are  of 
independent  interest  in  computer  vision,  as  they  provide  new  methods  for  surface 
reconstruction  that  overcome  many  of  the  same  problems  that  had  confronted  researchers 
in  oceanography,  namely  computational  complexity  and  the  ability  to  calculate  error 
statistics. 

Finally,  one  of  the  other  strengths  of  these  multiresolution  models  is  that  they  also 
allow  us  to  calculate  likelihood  functions  for  our  models,  pemitting  us  to  both  validate 
models  and  estimate  parameters.  Applications  of  this  capability  in  both  oceanography 
and  computer  vision  are  given  in  the  references  cited. 
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